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1 INTRODUCTION 



ABSTRACT 

The gamma-ray burst (GRB) / X-ray flash (XRF) events GRB 031203, discov- 
ered by INTEGRAL, and XRF 060218, discovered by Swift, represent two of only five 
GRB-SNe with optical spectroscopic confirmation of their SN components. Yet their 
observed high-energy properties offer a sharp contrast: While GRB 031203 was de- 
tected as a short 40-s burst with a spectrum peaking at -Epoak > 190 keV, XRF 060218 
was a Tgo ^ 2100-s long, smoothly evolving burst with peak energy E'pcak — 4.9 keV. 
At the same time, the properties of the two expanding dust-scattered X-ray halos 
observed in a fast-response XMM- Newton observation of GRB 031203 reveal that this 
event was accompanied by an "X-ray blast" with fluence comparable to or greater 
than that of the prompt gamma-ray event. Taking this observation as our starting 
point, we investigate the likely properties of the X-ray blast from GRB 031203 via 
detailed modeling of the XMM data, discovering a third halo due to scattering off 
a more distant dust sheet at ds = 9.94 ± 0.39 kpc, and constraining the timing of 
the X-ray blast relative to the GRB trigger time to be to — 11 zL 417s. Using our 
constraints, we compare the properties of GRB 031203 to those of other GRB-SNe 
in order to understand the likely nature of its X-ray blast, concluding that a bright 
X-ray flare, as in GRB 050502B, or shock breakout event, as in XRF 060218, provide 
the most likely explanations. In the latter case, we consider the added possibility that 
XRF 060218 may have manifested an episode of bright gamma-ray emission prior to 
the burst observed by Swift, in which case GRB 031203 and XRF 060218 would be 
"cosmic twin" explosions with nearly identical high-energy properties. 

Key words: dust, extinction — gamma-rays: bursts — supernovae: general — su- 
pcrnovae: individual: SN 20031w — supernovae: individual: SN 2006aj 

anomalous in their own way; even considered as a group, the 
properties of GRB-SNe thus resist ready generalisation. 



The discovery of the connection between long-duration 
gamma-ray bursts (GRBs) and X-ray flashes (XRFs), on 
the one hand, and Type Ibc supernovae (SNe) on t he other, 
has co nfirmed the collapsar mo del as proposed by WooslevI 



( 199i ) - for a recent review, see IWooslev fc BloomI j20oi ^ 



Several distinct types of investigation provide support 
for the GRB- (XRF-) SN connection; however, the most 
convincing evidence has come from spectroscopic observa- 
tions of coincident supernovae associated with the lowest- 
redshift GRBs and XRFs. As is frequently the case in astron- 
omy, these nearby and best-studied cases have each proven 



In all, five spectroscopically confirmed GRB-SNe - a 
term we will use to refer to the GRB- and XRF-SNe, as 
a group - have been observed to-date. SN 1998bw offered 
the first direct evidence for GRB-SN association, as the 
SN that appeared in the error box of GRB 980425, with 
an explosion time consistent with the GRB trigger, was 
shown to be excep tionally luminous and energetic in the 
optical (jGalama et al. 1998). with an associated mildly rel- 
ativistic (F >, 3) outfl ow that made it uniquely radio-bright 
( Kulkarni et al.iri998l \ Occurring in a catalogued galaxy at a 
distance of just 36Mpc {z = 0.0085), however, GRB 980425 
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was four orders of magnitude less luminous than the typical 
z > 1 burst. 

XRF 020903, at z = 0.25 the first XRF with a 
determined redshift and observed optica l and radio af- 



terglow emission (jSoderberg et al.l l2004al'). also provided 



the first XRF-associated supernova IjSoderberg et al.l 12005 



Bersier et al ; its spectru m was ultimately shown to be 

similar to that of SN 1998bw jSoderberg et ahlbood ). con- 
firming the deep connection between XRFs and GRBs. 

GRB 030329 / SN2003dh at z = 0.17 provided the first 
definitive association of a "cosmological" GRB with a coinci- 
dent type Ic supernova, since it had an is otropic-equivalent 
gamm a-ray output of iJ^.iso ~ 10^^ erg ( Vanderspek et al 
20041 ). within the observed range for z ^ 1 events. Broad- 



line S N Ic features were visible within a few days of the 
GRB (|stanek et al.ll2003l ). and the ultimate evolution of 
the supernova showed a striking similarit y to SN 1998bw 



I Hiorth et al 



20031 : 



Matheson et al 



20031 ). 



two earlier associations. 

Later that year, GRB 031203 was detected by IN- 
TEGRAL with a duration of 40-s and a peak e nergy 



of > 190 keV (ISazonov et all |2004 hereafter [SLSOJ). Six 
hours after the burst, a fast-response XMM-Newton ob- 
servation of the burst position began; a fading X-ray 



after glow was quickly identified (jRodriguez-Pascual et al 
2003), coincident with a, z = 0.105 galaxy (jProchaska et al 
200J). This was confirmed as the host galaxy once an 



associated type Ic s upernova, SN 20031w, was discovered 
( Malesani et al.ll2004l ). Investigation of the radio and X-ray 
energetics of the event revealed further similarities to 
SN1998bw dSoderberg et al.ll2004bl ). althou gh the SN itself 



was s omewhat subluminous by comparison (|Gal-Yam et al 
20041) . 



Contemporaneously, analysis of the XMM data revealed 
two exp anding X-ray halos ce ntred on t he fa ding GRB af- 
terglow (|Vaughan et al.l |200J, hereafter |V0J) ; these halos 
were convincingly interpreted as scattered photons from 
a bright X-ray "blast" that occurred near-simultaneously 
with GRB 031203. Mode lin g the pro perties of the X-ray 



bl ast.lWatson et al.1 l|2004l . or lW04h and lWatson et al.1 liooi 



WOd ) argued that the prompt emission of GRB 031203 



was, in fact, an X-ray flash, exhibiting greater fluence in 
X-ray compared to gamma-ray, log(g 2-30kev/»5'30-400kev) = 
0.6±0 .3 > 0. Ind ependent analysis bv lTiengo fc Mereghetti 
( 20061 . or 'tMOG') yielded a lower X-ray fluence estimate 
which nonetheless also satisfies S2-30keV > 53o-4ookoV. 

The nature of GRB 031203 is not easily resolved, for 
two reasons. First, the bright soft X-ray emission cannot 
be interpreted as the low-energy tail of GRB 031203, as its 
fluence significantly exceeds t he low-energy extrapo lation 
of the INTEGRAL spectrum (|SLS04l : IW06l : iTMod ). Sec- 
ond, the INTEGRAL observations themselves would have 
detected the X-ray blast if it was emitted over a brief time 
interval prior to or within 300 s following the burst itself, be- 
fore 22:06:27 UT, when INTEGRAL began to slew to a new 
pointing. Thus, two distinct explosive episodes are required, 
with the X-ray blast most likely occurring both physically 



and tem porall y apart from the GRB. To satisfy these con- 
straints, Iwod consider the possibility that GRB 031203 was 
followed by a bright X-ra y flare such as that observed by 
Swift from GRB050502B (jBurrows et al.ll2005h : such fiares 
- albeit rarely as bright as that case - are familiar features of 
the GRB X-ray lightcurves now routinely gathered by Swift 



(e.g., IChincarini et al.ll2007l ') 



XRF 060218, detected by Swift (ICampana et al 



2006) 



and coincident with a z = 0.033 host galaxy ([Mirabal et al 
2006[ ). evolved into a type Ic supernova that was ob- 



served from its earlies t moments, including extensive spec- 
troscopic observatio ns (jMirabal et al ]|2006l :l iPian et al ]|2006l : 



Mazzali et al.l 20061 ). thanks to the prompt alert and ex- 
tended follow-up observations of Swift. Analysis of the 
X-ray, optical, and radio energetics of XRF 060218 demon- 
strate a deep similari ty to SN 1998bw and GRB 031203 



vahdatmg the jSoderberg et ahlbood ): however, XRF 060218 remains an 



anomaly in that its prompt emission exhibited a uniquely 
extended evolution, having a 90%-containment duration of 
Tgo — 2100 s and a correspondingly low peak fiux Fpcak = 
8.8x10"^ erg cm~^ s'^ (15-150 keV). This peak fiux is more 
than an order of magnitude lower than those of all other 
GRB-SNe except for XRF 020903, which had an even softer 
spectrum; indeed, it was only detected because of the "im- 
age trigger" capability of the Swift BAT, which identified 
XRF 060218 as an uncatalogued source apparent in the first 
64 s integration taken at that pointing. 

Emission from XRF 060218 is present in the first image 
taken by Swift on that particular orbit; as such, a GRB sim- 
ilar to GRB 031203 could have occurred prior to the start 
of Swift observations. Moreover, a significant fraction of the 
prompt X-ray emission of XRF 060218 is supplied by a ther- 
mal component, unobserved in other GRB-SNe, that is at- 
tributed to the shock breakout of the associated supernova 
( Campana et al.ll2006l ) rather than to the relativistic jet usu- 
ally invoked to explain the prompt emission. If the X-ray 
blast from GRB 031203 can be explained as a similar shock 
breakout, then these two events - despite their very differ- 
ent appearances - would in fact be "cosmic twins" exhibit- 
ing nearly iden tical evolution in the X-ray and gamma-ray 
bands (see also Ghisellini et al. 20061 ). 



In this paper, we undertake a quantitative exploration 
of this hypothesis by constraining the properties of the X-ray 
blast in GRB 031203 via detailed modeling of its expanding 
X-ray halos; this analysis also allows us to address the al- 
ternate theory that the X-ra y emis sion from GRB 031203 
was due to a later X-ray flare IWOel ). §2 presents the XMM 
observations of the afterglow and halos of this burst and 
the details of our analysis, including derivation of param- 
eter uncertainties. §3 discusses our results in light of vari- 
ous hypotheses regarding GRB 031203 and XRF 060218, and 
places these events in the larger context of the known vari- 
eties of cosmic explosion. §4 summarizes our conclusions. 
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2 OBSERVATION AND DATA ANALYSIS 

The GRB 031203 data (Observation ID 0158160201), down- 
loaded from the XMM-Newton archive, are processed and 
reduced for EPIC MOS and PN cameras through the stan- 
dard procedures (tasks emchain and epchain), with stan- 
dard filtering, using the XMM Science Analysis System 
(SAS), version 7.1.0. The observation had start and stop 
times of 04:10:18 UT and 20:16:20 UT on 4 Dec 2003, which 
after standard filtering yields 57.8 ksec of good time and 
57.3 ksec net exposure in CCD 1 of EPIC-MOSl and M0S2, 
and 56.2 ksec of good time and a CCD-dependent net expo- 
sure between 53.4 and 54.0 ksec in EPIC-PN. The observa- 
tion was not affected by significant soft-proton fiaring. 



to increasing Dec, as usual. The temporal coordinate t is 
measured relative to the burst trigger time (t = 0), scaled 
such that the XMM observation begins ai t — 1, 22.2 ksec 
(6'^10™44'^) after the burst. Using this temporal coordinate, 
the XMM observation extends from t = 1 (defined by the 
earliest arrival time in our dataset, from the MOS-2 detec- 
tor) to t m 3.61. An illustration of the resulting dataset 
is presented as an image in normalized coordinates {x/y/i, 
y/Vt) in Fig.IU 

We represent the fading GRB afterglow using an ana- 
lytical model for the XMM Point Spread Function (PSF), 
which approximates the PSF as a King function: 



(l + rVr?)- 



(1) 



2.1 Model Construction and Fitting 

We consider a circular region with radius 311" centred on 
the GRB afterglow, and restrict our attention to detected 
counts in the energy range 0.2 to 3.0 keV (0.2 to 4.0 keV 
for PN). We generate an exposure map with 1" pixel res- 
olution, using the eexpmap task, to account for PN active 
regions; exposure map values are treated as binary (on or 
off) in the subsequent analysis, and no initial exposure map 
is created for the MOS data as our circular region of interest 
is fully contained on the central CCD. Bright X-ray sources 
within the region of interest are identified and all counts 
within an exclusion radius of 15" of each source position are 
excluded from analysis; these exclusion regions are then in- 
tegrated into exposure maps for the PN and MOS data for 
analysis purposes. After all exclusions we accept for analysis 
10,754 counts from the MOS-1 detector, 10,535 counts from 
MOS-2, and 32,085 counts from PN. 

We construct a 16-parameter, three-dimensional (x, 
y, t) maximum-likelihood model including parameters to 
describe the fading GRB afterglow, the expanding dust- 
scattered halos, and the vignetted X-ray (and charged 
particle-induced) background. The total counts from all 
model components, for the 2-MOS and PN datasets, are 
constrained to be equal to the observed count totals given 
above. The model probability density is evaluated at the 
position and arrival time of each detected photon, the value 
of the probability density function is interpreted as being 
proportional to the likelihood of observing that event, and 
a global likelihood calculated; we then attempt to maximize 
the global likelihood using numeric al techniques, spec ifically, 



the AMOEBA minimization routine (Press et al 



19951 ) 



plemented in the IDL programming environment. For ease 
of modeling, we convert photon angular positions and arrival 
times as follows: angular coordinates x and y are measured 
in arcsec relative to an initial estimate of the GRB posi- 
tion, XMM physical coordinates (26024.8, 23674.1), with 
the +x direction corresponding to decreasing R.A. and +y 



where a is the King slope, is the King core radius, r is the 
radial distance from the source (ring centre), and the inte- 
grated surface brightness is normalized to unity for a > 1. 
Since we combine the data from the two MOS cameras, we 
use average MOS PSF parameter values for the MOS data 
(r^ = 4.916 and a = 1.442), and PN PSF values for the PN 
datE0 (re = 6.636 and a = 1.525). The fading of the after- 
glow is parameterized as a power-law with temporal index 
ag, referenced to the burst trigger time. 

We represent the expanding, fading dust-scattered ha- 
los, centred on the GRB position, using a numerical model 
for ring-like sources observed with XMM that we have de- 
veloped. The derivation of this model is provided in Ap- 
pendix [X] as realized in our code, this bivariate function 
is precalculated on a grid and interpolated to the particu- 
lar values required for each function evaluation. Fading of 
the halos is parameterized as a power-law with temporal 
index a^, referenced to the burst trigger time; this index 
is assumed to be the same for all rings. In order to avoid 
complications associated with normalizing the model prob- 
ability density function for indices < — 1 as well as for 
Oh > —1, we impose a hard limit ah > —0.95 on this param- 
eter during fitting. Initial investigations suggested that this 
limit was encountered in less than 5% of bootstrap trials; 
however, our final analysis resulted in a somewhat greater 
proportion of 9% of all trials encountering this limit. 

Because the rings expand by a factor of \/3.6 ~ 1.9 over 
the course of the observation, and because we wish to ac- 
curately model the background over the region of interest, 
we find it necessary to apply an off-axis vignetting function 
i Lumb et al.. ,2003') to the rings and background: 



5(r) = S(0)(l- 



(2) 



where 5(0) is the surface brightness on-axis, S{r) is the sur- 
face brightness at off-axis distance r (measured in arcsec), 
and m and e are the vignetting parameters for the EPIC 
cameras. We use m = 3.43 X lO"'^ and e = 1.45 for the MOS 
cameras, and m — 3.51 x 10~^ and e — 1.53 for the PN, 



^ PSF model and parameter values are obtained from the XMM calibration document "In Flight Calibration of the PSF for the PN 
Camera" by Simona Ghizzardi, located at http://xmm.vilspa.esa.es/docs/documents/CAL-TN-0029-l-0.ps.gz See, in particular, Eq. 1 
and Table 1 in the document. 
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Figure 1. (left) Images of the 2-MOS and PN counts data in normalized coordinates {x/\/t, yj^Jt). The fading GRB afterglow, 
Halo 1 (outer), and Halo 2 (inner) are readily distinguished; the signature of the more distant and fainter Halo 3 (close to 
the GRB) is more subtle. White streaks in the images are due to bright source exclusion regions and PN camera chip gaps, 
which are incorporated into our model via MOS and PN exposure maps. 



Figure 2. (right) Histogram of 2-MOS + PN counts data (black), best-fit model (green), and residuals (lower panel), binned 
according to normalized radius, r jyfi. Contributions of the fading GRB afterglow (GRB) and the three expanding halos (HI, 
H2, H3) are individually distinguishable; the fall-off in counts at large normalized radius is due to our region of interest cut 
(''max = 311") as smeared out by the normalization procedure. 



which we find accurately models the vignetting of PN and 
MOS exposure maps generated in SAS. 

The excluded sources and the PN exposure map present 
discontinuities to the model probability density function, 
which must be normalized to yield an appropriate fit. We 
perform this normalization via numerical integration over a 
three-dimensional rectangular grid (x, y, t), where the grid 
spacing is 2" in x and y, and the temporal grid is divided 
into 40 equal intervals. The same cuts applied to the data, 
including the circular region of interest, the exclusion regions 
around identified bright sources, and (for PN only) the PN 
exposure map, are applied to the grid prior to integration, 
and are thus (approximately) accounted for. 

We maximize the likelihood function using the AMOEBA 
algorithm, in stages, using multiple fitting rounds at each 
stage (with fit starting points randomized prior to each 
round) to ensure a relatively wide-ranging exploration of the 
parameter space. In sequence, we execute: (1) Ten rounds of 
fitting the GRB counts parameters and fading index Og-, (2) 
Ten rounds of fitting with all ring parameters (counts, radii, 
and fading index q^) also free; and finally, (3) Five rounds 
of fitting with all parameters free. 

Initial examination of the data suggested the presence 
of a third dust-scattered halo at small radius (Fig. [l]). Pre- 
liminary investigation confirmed the reality of this feature, 
and indicated that it should be interpreted as due to X-ray 
scattering off a third dust sheet, more distant than the two 
discussed by previous authors ( i]3.2|) . We therefore incorpo- 
rate this third halo into our fits, with the same burst time 
to and fading index as the others. 

Uncertainties in our model parameters are derived with 
the bootstrap Monte Carlo method, executing multiple trials 
and accumulating best-fit parameter values from each trial 
to realize the posterior distribution of our model parameters. 



For each trial, a new dataset is generated by drawing ran- 
domly, with replacement, from the true dataset until we have 
the same number of events as in the actual data. The fitting 
sequence, described above, is identical for each trial; param- 
eter starting points are randomized based on preliminary 
estimates (subsequently verified in our final analysis) of the 
uncertainties in each parameter. The result should approx- 
imate the posterior distributions that would result from a 
full Bayesian analysis using uniform (non-informative) prior 
distributions for each of the model parameters. For our final 
analysis, presented below, we executed 1000 trials in all. 

The optimization of nonlinear functions in many- 
dimensional spaces is a known hard ("NP") computational 
problem. We therefore point out that the bootstrap Monte 
Carlo approach is expected to yield conservative estimates 
of parameter uncertainties in cases where the model fails to 
offer a statistically-complete representation of the data, or 
where the underlying fitting procedure fails to identify the 
global minimum in each and every trial. The reason is that, 
quite simply, the failure to identify the best possible fit is 
robustly expected to drive the fitting routine to explore a 
larger (rather than smaller) region of parameter space. This 
is true as long as the fitting routine is not trapped in a single 
local minimum, repeatedly, in trial after trial; in our analysis 
we avoid this numerical disorder by randomizing parameter 
starting values for each trial, as mentioned above. 

2.2 Results 

The results of our model fitting, including parameter un- 
certainties quoted both as standard deviations and as 
90%-confidence intervals, are presented in Table [U 90%- 
confidence intervals are defined as the minimum-length in- 
tervals that encompass 90% of the parameter values from our 
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bootstrap Monte Carlo trials (e.g., 900 of 1000 parameter 
values). Associated background count totals can be derived 
by subtracting the counts of the model components given 
in the table from the total instrument counts of 21,289 for 
2-MOS and 32,085 for PN (SgHJ. 

To illustrate the model goodness-of-fit, we present in 
Fig. [2] a histogram of the data (2-MOS + PN counts) com- 
pared to our best-fit model. In this figure, as well as in Fig. [1] 
we present the data in terms of the "normalized radius," 
which is the radial distance from the GRB afterglow po- 
sition, divided by the square root of the normalized time 
(time measured from the GRB trigger, with t = 1 at the 
start of the XMM observation). This normalization serves to 
highlight the expanding halos as distinct features, while not 
significantly altering the appearance of the afterglow itself 
(at small radius); the fall-off in counts at large radius is due 
to our region of interest cut (rmax = 311"), as smeared out 
by the conversion to normalized units. We find the overall 
fit to be satisfactory. Specifically, although several system- 
atic features are apparent in the residuals, we have investi- 
gated these without identifying any obvious model failings 
that might account for them. In particular, we have veri- 
fied that none of the features are due to unaccounted-for 
bright point sources or instrument artifacts (e.g., hot pixels 
or hot columns). Variations in effective exposure within our 
region of interest, not accounted for by our simple vignetting 
model, may be responsible for some of these trends. Note 
that our fitting procedure is executed against the unbinned 
data, in three dimensions, and that this figure merely rep- 
resents one possible projection of the model and data onto 
a one-dimensional space, for visualization purposes. 

In Fig. [3] we present the full posterior distribution for 
four parameters of special interest: Og, the temporal power- 
law index for the fading of the GRB afterglow; to, the time of 
zero-radius for the expanding halos; ah , the temporal power- 
law index for the fading of the expanding halos; and the sum 
of the 2-MOS and PN counts for Halo 3, the halo due to the 
most distant scattering dust sheet, discovered in our anal- 
ysis. 90%-confidence intervals on each of these parameters 
are indicated. 



3 DISCUSSION 

In the sections below we discuss the prior understanding of 
the GRB 031203 X-ray blast phenomenon and our own re- 
sults as to the properties and timing of the event. We then 
discuss the full set of observations of GRB 031203 in the 
context of the known varieties of cosmic explosion, includ- 
ing XRF 060218. 



3.1 Previous Findings 

Ivoi originally reported the discovery of two evolving dust- 
scattered X-ray halos centred on GRB 031203 and made 
the first attempt to model the halos and extract parame- 
ters of the X-ray blast and associated dust scattering sheets 



along the line of sight within the Milky Way. They con- 
strained the time of burst for the inner and outer rings to 
be to = 2794:tll^l s and to = 2005tlll7 a respectively, rela- 
tive to the INTEGRAL trigger, and by applying a model of 
the dust properties, determined that the initial soft X-ray 
event had a power law spectru m wit h phot on index F ~ 2. 

Subsequently, |W04|. Iwod . and iTMod . using different 
analysis approaches, refined the details of this picture. They 
derived more precise distances to the dust sheets (di « 
870 pc and d2 ~ 1390 pc, respectively), and correspond- 
ingly more precise constraints on the timing of the X-ray 
blast. The time constraint from WOd is to = 600 ± 700 s at 
90% confidence; the inferred power-law spectral index for 
the X-ray blast is repo rted as either F — 2.5 ± 0.3 (|W0J) or 
F = 2.1 ± 0.2 iTMod l. results which are consistent at the 
~la level. The two conclusions as to the fluence of the X-ra: 

2 



blast over 1-2 keV, Fx = 1.5(31 x 10 
versu s 3.6(2) x 10~^ erg cm" 



' erg cm (jWOi 



7 -2 jTMod ). are more discrepant 



which lTMOg attribute to the authors' diverging assumptions 
on the relation between X-ray scattering optical depth and 
extinction. 

In discussing their model for the dust properties, Wod 
show a plot of the halo fading (their Fig. 2) which cor- 
responds to a temporal power-law index of ~ —0.6. 
They use this fact (along with the size of the scatter- 
ing halos) to derive a maximum size for the dust grains, 
ctmax = 0.50±0.03 ^m; in accounting for the effects of smaller 
grains they use a power-law size distribution with power-law 
index equal to —3.5. 

3.2 Discovery of a Third Halo 

Our analysis reveals the presence of a third, faint halo 
(Figs. [1] [2)| in addition to the two identified by previous 
authors. The presence of this halo is robust in our analy- 
sis, with associated total 2-MOS -I- PN counts between 403 
and 572 at 90%-confidence (Fig.|3}i), and with roughly equal 
counts in the two sets of detectors, which are fit indepen- 
dently (Table [TJ. 

Initially, we speculated that the third halo might be 
the signal of a second X-ray blast or early afterglow emis- 
sion, with the X-rays scattering off one of the two previously 
known dust sheets. If this were the case, the third halo would 
necessarily be associated with a distinct (later) to and would 
share the same distance d as one of the other halos. 

To investigate this possibility, we performed a distinct 
bootstrap Monte Carlo analysis of 600 trials, allowing to for 
the third halo to vary while fixing to — for the two other 
halos. Given the angular size of the halo at the start of the 
XMM observation, 9, to is directly related to the distance to 
the dust scattering sheet by d = 2c{t — to)/ 6'^, where t — to 
is the time from the burst to the start of XMM observa- 
tions. This analysis yields a 90%-confidence interval for ds 
of 8597 to 10503 pc, inconsistent with both di « 870 pc and 
d2 ~ 1390 pc. Hence, we conclude that the third halo is due 
to a third dust sheet and require to to be identical, for all 
halos, in our full analysis. 
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Table 1. Fit parameters for GRB 031203 afterglow and expanding halos 



Parameter 


Median 


cr 


90%-conf. 


GRB X 


0.46 


0.14 




U. / U 


y 


0.19 


0.13 


— U.U4 


n /I n 
U.4U 


2-MOS counts 


2252.7 


60.3 


ZiO I .o 




PN counts 


2340.0 


125.5 




ZoOy.Z 


Fading index, cxg 


-0.557 


0.056 


— U.u4 


— U.4D 


Halos to (s) 


11.0 


416.9 


— DzU.D 


i ZZ.o 


Fading index, 


-0.77 


0.10 


1^ — u.yo ) 


— U.Do 


Halo 1 Radius (arcsec) 


145.04 


0.70 




14D.U ( 


Distance (pc) 


871.7 


9.7 


855.6 


886.8 


2-MOS Counts 


1365.0 


67.8 


1251.1 


1467.4 


PN Counts 


1686.4 


88.8 


1568.4 


1855.6 


Halo 2 Radius (arcsec) 


114.93 


0.63 


113.90 


115.94 


Distance (pc) 


1388. 


14. 


1364. 


1410. 


2-MOS Counts 


2059.3 


65.5 


1952.9 


2164.1 


PN Counts 


2229.3 


108.9 


2041.6 


2389.5 


Halo 3 Radius (arcsec) 


42.97 


0.81 


41.46 


44.06 


Distance (kpc) 


9.94 


0.39 


9.32 


10.52 


2-MOS Counts 


268.5 


40.7 


203.4 


334.1 


PN Counts 


219.7 


32.2 


177.7 


274.8 



Median, standard deviation about the median (a), and 90%-confidence 
intervals are derived from model fits to 1000 bootstrap Monte Carlo re- 
alizations of the dataset. PN and 2-MOS counts arc defined as the total 
counts gathered from the PN (0.2 to 4.0 keV) and two MOS cameras (0.2 
to 3.0 keV), respectively. Halo radii are defined as the radius in arcsec 
at the start of XMM observations. Halo distances are determined by the 
best-fit halo radii and ip values from each trial and are not independent 
parameters. GRB coordinates (x, y) are offsets in arcsec from XMM phys- 
ical coordinates (26024.8, 23674.1); the median offset of (0.46, 0.19) arcsec 
corresponds to R.A. 08'^02™3W20, Dec. -39°51'02?8 (J2000) using the 
uncorrected XMM aspect, to is the time of zero-radius for the expand- 
ing halos relative to the INTEGRAL trigger time of 22:01:28 UT. The 
model enforces a hard limit > —0.95 which is encountered in 9% of our 
bootstrap trials. 



The third dust sheet is located at greater distance and 
greater Galactic radius [R ~ 14.5 kpc) than the two nearby 
sheets, and nominally at greater distance from the Galac- 
tic plane as well {z ~ —0.83 kpc), given the Galactic co- 
ordinates toward GRB 031203, I = 255.7° and b = -4.8°. 
However, reconstructions of the three-dimensional distri- 
butio n of neutral hydrogen density throughout t he Milky 
Way jKalberla et alJboOTl : iKalberla fc Dedeslboosi ) reveal a 
"warp" in the disk which has an amplitude of 5^ « 1 kpc at 
R = 15 kpc, in precisely the {—z) sense needed to place the 
third dust sheet within or near the densest portion of the 
Hi disk. The existence and location of the third dust sheet 
is not particularly surprising in this context. 

Our final analysis yields 90%-confidence intervals on 
the distances to the dust sheets as follows: d\ — 856- 
887 pc, d2 = 1364-1410 pc, and = 9321-10523 pc (Ta- 
ble [T|; these uncertainties are inclusive of uncertainties in 
the value of to itself, which we have not fixed (e.g., Fig.[3}3). 
For Halo 1 and Halo 2, our values agree with those derived 
by Iwod (di = 868ligpc and da = 1395+30 pc) and ItM06| 
(di = 870 ± 5pc and da = 1384 ± 9pc); note, however, that 
TMod fix to = s in their analysis. 

The median value of the halo fading index from our 
analysis is ~ —0.77 with a 90%-confidence interval 



Uh < —0.65 that is unconstrained from below because of 
our hard limit {an > —0.95); Fig. |3j; presents the full poste- 
rior distribution. Our result appe ars mildly in confiict with 
the value an ~ —0.6 from WOfil . with associated implica- 
tion s for th e properties of the scattering dust particles (see 
also TMod ): we have not investigated these implications in 
detail. 



3.3 Properties of the GRB 031203 X-ray Blast 

Our analysis yields a 90%-confidence interval on the timing 
of the X-ray blast which places it within roughly ten min- 
utes of the GRB 031203 trigger, between —621s and +723 s, 
with median value and standard deviation to — 11 i: 417 s; 
the full posterior distribution for to is shown in Fig. ^jp. 

Our constraint is comparable in precision to the best 
previo us determination, to = 600 ± 700 s at 90%-confidence 
( W06l ). The shift of our confidence interval to earlier times 
sharpens the contrast between the implications of the XMM 
and INTEGRAL observations; The INTEGRAL pointing 
extends to 300 s after the GRB trigger, and only 22.5% of 
our trials yield to values consistent with this constraint. If we 
assume, in accordance with standard GRB models, that the 
X-ray blast could not have preceded the gamma-ray trigger, 
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Figure 3. (left) Approximate posterior distributions of selected model parameters, as derived from our bootstrap Monte Carlo 
analysis, (a) The GRB afterglow temporal fading index, oig; (b) to for the expanding halos, corresponding to the characteristic 
emission time of the X-ray blast, measured relative to the trigger time for GRB 031203; (c) The halos temporal fading index, 
ah (note this parameter is subject to a hard limit, oth > —0.95); (d) The sum of 2-MOS + PN counts for the expanding 
Halo 3, discovered in our analysis. 90%-confidence intervals for each parameter are indicated graphically (dotted lines). 



Figure 4. (right) Superposition of XRF 060218 Swift BAT (green, 15-350 keV) and XRT (red, 0.3-10 keV) light curves with 
a hypothetical precursor event having properties identical to GRB 031203 (black, 2It2QQ,keV). The start of Swift BAT obser- 
vations defines t = 0; the INTEGRAL light curve of GRB 031203 is obtained from[SLSOJ, rebinned, corrected for redshift and 
time dilation, and offset by —70s to represent its occurrence prior to the start of Swift observations. Negative (background- 
subtracted) flux measurements are not shown on this logarithmic plot. 



then our upper limit of to < 545 s at 90%-confidence is also 
relevant in this context. 

An alternative approach to resolving this dilemma 
would reduce the peak flux of the X-ray blast - making it 
undetectable in the INTEGRAL observation - by extending 
its fiuence over a significant time interval after the GRB. 
This hypothesis is also constrained by our analysis, since 
we observe no evidence that the halos are radially resolved; 
our analytical model for unresolved ring-like sources (Ap- 
pendix [A]) provides a satisfactory fit to the radial profiles of 
the halos (Fig.EJ. 

While we have not constrained the duration of the X-ray 
blast directly, our results on the timing of the blast demon- 
strate that our model is sensitive to changes of Sto ~ 600 s, 
and we expect that, in demonstrating an adequate fit with 
the unresolved "ring" model, we have constrained the du- 
ration of the blast at approximately this magnitude, St <. 
600 s. Given that the XMM observation began t = 22.2 ksec 
after the GRB trigger, a duration of 600 s corresponds to 
5t/t — 3% or Sri ~ 4" for Halo 1, which would be most 
sensitive to these effects. Since the XMM PSF has a width 
of FWHM — 4", it is reasonable to think the data would be 
sensitive to signs of broadening at this level. Moreover, the 
dust sheets themselves must be extended to some degree; 
any such extension would tend to further broaden the radial 
profile of the halos, beyond any broadening due to the finite 
duration of the blast itself. 



Finally, we note that |W04| call attention to the relatively 
slow fading of the X-ray afterglow of GRB 031203 during the 
XMM observation, ag ~ —0.55 (Table [T]), which was signif- 



icantly flatter than the canonical ax ~ —1.3 familiar from 
Beppo-SAX observations at the time. In retrospect we can 
identify this as one of the first observations of a "plateau 
phase" in GRB X -ray afterglows, now familiar from Swif t 
observations (e.g., Nousek et al. 20061 : Racusin et al. 20091 ) . 
As we have since learned, extrapolation of this slow decay 
to early times, t <; 1000 s, is not particularly revealing of the 
likely X-ray flux at that time. 



3.4 Could GRB 031203 and XRF 060218 be 
Cosmic Twins? 

As mentioned in iJT] a GRB similar to GRB 031203, as ob- 
served by INTEGRAL, could have occurred immediately 
prior to Swift observations of XRF 060218 without being 
observed by Swift. Likewise, a slow-evolving XRF similar 
to XRF 060218 could have occurred following GRB 031203 
without being observed by INTEGRAL, if it exhibited either 
a sufficiently low peak flux or peaked after t = 300 s, when 
INTEGRAL began its slew to a new pointing. In this sce- 
nario, then, both events would consist of two separate bursts, 
one peaking at gamma-ray energies and one peaking in 
the X-ray. Importantly, while we have evidence for distinct 
gamma-ray and X-ray emission episodes in GRB 031203, 
in this picture only the X-ray episode was observed from 
XRF 060218. 

A variation on this scenario, without any separate 
episode o f gamma-ray em i ssion from XRF 060218, was pro- 
posed bv iGhisellini et all (|2006l '). These authors were con- 
cerned with the question of whether the prompt high-energy 
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emissions of GRB 031203 and GRB 980425 could be recon- 
ciled with the high-energy relationships between peak energy 
and isotropic-equivalent energy output that are observed 
among the z >1 GRB population. They concluded, as we do 
here, that the observation of XRF 060218 by Swift suggests 
strongly that a similar phenomenon might be responsible for 
the "X-ray blast" inferred from the expanding X-ray halos 
of GRB 031203. 

To illustrate the nature of our proposal, we have gen- 
erated a figure that su perimp oses the prompt-emission light 
curve of GRB 031203 (jSLSOJ), appropriately rescaled to the 
redshift and luminosity distance of XRF 060218, on the Swift 
BAT and XRT lightcurves of the prompt and early emission 
of XRF 060218 (Fig. g]). With t = at the beginning of 
Swift observations, it is necessary to offset the INTEGRAL 
light curve to earlier times in order to avoid detection by 
the Swift BAT; we arbitrarily choose a —70s offset. The 
BAT light curve is a combination of the mask- weighted light 
curve generate d by the BAT standar d analysis {t < 350 s) 
with data from Campana et al. ( 20061 ) (64-s integrations at 
t > 350s). The XRT light curve (red) for XRF 060218 is 
obtained from the Swift/ 'KRT light curves reposi tory at the 
University of LeicesteJl (|Evans et al.ll2007l . \200d ). 

As Fig. g] shows, this "cosmic twin" to GRB 031203 
would have exhibited a peak flux of 2.7x 10~^ erg cm~^ s~^ 
(20-200 keV), roughly 11 times brighter than the origi- 
nal event, which occurred at about three times the dis- 
tance of XRF 060218. This addresses an important con- 
dition on our hypothesis, which is that the burst not be 
bright enough to t rigger detectors of the Interplanetary 
Network flPN: e.g.. iHurlev et allbood '). and in particular. 



Konus-WIND (|Aptekar et aLlUQm which did not report 



any trigger associated with XRF 060218. A review of the 
catalog of Konus-WIND short burst detection^ reveals that 
the peak flux distribution peaks at 7.5xl0~® erg cm~^ s~^ 
(15-2000 keV), although bursts are detected down to a min- 
imum peak flux of 2x10^^ erg cm~^ s~^. Depending on the 
high-energy spectrum of GRB 031203, the peak flux of its 
hypothetical z — 0.033 cosmic twin would be from 5.2x 10~® 
to 7.9x10"" erg cm"^ s"^ (15-2000 keV), likely below the 
regime of maximum sensitivity for this experiment. With 
a peak flux falling close to threshold in the brightest case, 
relatively slight variations in the profile or spectrum of the 
event would be sufficient to render the burst undetectable. 

An associated "classical GRB" like GRB 031203 thus 
seems plausible; however, our analysis indicates that the 
X-ray properties of the two events must differ to some 
degree. First, we have constrained the timing of the 
GRB 031203 X-ray blast to fall within to = H ± 417 s of 
the GRB (to < 545s at 90%-confidence; ^3.3p . whereas the 
XRT lightcurve for XRF 060218 rises slowly until it peaks at 
t ~ 1000 s. Second, depending on the assumed rela tion be - 
tween extinction and X-ray scattering optical depth ( TM06t ). 
the GRB 031203 X-ray blast may have released as much 



as three times the energy of the prompt X-ray emission of 
XRF 060218 (Table[2|. 

To illustrate these constraints in detail, and explore 
other possible explanations for the GRB 031203 X-ray blast. 
Table [2] presents the gamma-ray and X-ray properties of all 
five confirmed spectroscopic GRB-SNe, and the properties 
of two other cosmic explosions that may be interesting in 
this context, GRB050502B and SN2008D. 

First, we consider the "cosmic twins" hypothesis, com- 
paring the properties of GRB 031203 and XRF 060218. 
A distinct, earlier episode of high-energy emission for 
XRF 060218 is readily accommodated by the total 
gamma-ray energy budget of GRB 031203; in order to avoid 
detection by various IPN experiments, the peak luminos- 
ity, peak energy, or total energy release should be some- 
what reduced from GRB 031203 values. The X-ray blast 
from GRB 031203 released more energy in the X-ray than 
XRF 060218; however, depending on the assumed relation 
between extinct ion and X-ray sc atterin g optical depth, the 
additional 25% IflMod ) to 200% jWOet ) energy requirement 
is not necessarily severe. Finally, the time delay from burst 
trigger to X-ray peak (technically, to the characteristic time 
of emission) should be at least twice as long for XRF 060218 
as for GRB 031203. None of these requirements appear to us 
prohibitive; collectively, however, they restrict the allowed 
parameter range of models considerably. The main appeal 
to satisfying these constraints, in an Occam's Razor sense, is 
the requirement for one fewer sui generis type of GRB-SN. 

Next, we consider th e X-r ay fiare explanation for the 
GRB 031203 X-ray blast (jwol). Even at the maximum in- 
ferred energy of 4x10''^ erg (2-10 keV), the ratio of X-ray 
to gamma-ray energies for GRB 031203 would be no greater 
than observed in GRB050502B. Moreover, the relative tim- 
ing of the gamma-ray and X-ray episodes is very close to 
what is needed to satisfy both our time constraint, derived 
from the expanding halos, and the INTEGRAL slew time 
constraint, t > 300 s. The latter applies to any relatively 
short (hence high peak fiux) model for the X-ray blast. At 
the same time, the X-ray flare from GRB050502B remain s 
an extreme example of its class fe.g.. lChincarini et al ] l2007l l. 
and in this model the GRB 031203 X-ray blast would be a 
similarly extreme case; in addition, it would be the only 
known example of an X-ray flare from a GRB-SN. 

None of the other cosmic explosions listed in Table [2] 
are known to exhibit strong distinct episodes of gamma-ray 
and X-ray emission. Their properties are presented to illus- 
trate the diversity of the GRB-SN phenomenon, and the 
gamma-ray and X-ray luminosities that may generally be 
expected from these events on the basis of past experience. 
Thus, a prompt high-energy emission episode in XRF 060218 
like those seen from SN 1998bw or XRF 020903, and occur- 
ring prior to the start of Swift observations, would certainly 
have fallen below IPN threshold, and if present, would help 
to homogenize the GRB-SNe as a class. 



^ Suij/t/XRT light curves repository : |http : //www ■ swift ■ ac ■ uk/xrt_curv6s/ 1 

^ Konus-WIND short burst catalog: [http : //www ■ lof f e ■ ru/LEA/shortGRBs/CatalogT] 
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Table 2. Properties of GRB 031203 and related cosmic explosions 



Event 2 logB-yiao log-Ex,iso Interpretation Refs. 



SN 


1998bw 


0.0085 


47.8 


47.1 


+10 


GRB-SN Ic 


1,2 


XRF 


020903 


0.251 


48.5 


48.9 





XRF-SN Ic 


3,4 


GRB 


030329 


0.167 


52.0 


51.1 


+3 


GRB-SN Ic 


5 


GRB 


031203 


0.105 


50.0 


49.2-49.6 


±600 


GRB-SN Ic + ? 


6,7,8,9 


GRB 


050502B 


(2.0) 


51.4 


51.0 


+250 


GRB + X-ray flare 


10,11 


XRF 


060218 


0.033 


49.4 


49.1 


+1050 


XRF-SN Ic + shock 


12 


SN 


2008D 


0.0065 


<46.5 


46.3 




SNIc shock 


13 



logE-y^iao is the isotropic-equivalent gamma-ray energy over 20 keV to 2 MeV, rest-frame (ergs); 
log i?x iso is the isotropic-equivalent prompt X-ray energy over 2 to 10 ke V, rest-frame (ergs) ; 



(s). A redshift of 2 = 2 is assumed for GRB050502B; no redshift-related corrections are ap- 
plied for S N 1998 bw and SN2008D. The two values of log^x.iao for GRB 031203 are from lTMod 
(49.2) and |W04| (49.6), respectively. The X-ray energy and time delay for XRF 06 0218 are de- 
rived by an alysis of the XRT lightcurve from the Swift XRT lightcurvcs repository ijEvans et al.l 
I2OO7I . I2OO9I ). The upper hmit on E^-^^o for SN2008D is extrapolated from the Sw ift BAT upper 
limit usi ng an a = 1.5 power-law spectrum. Refere n ces: fca lama et al.l lll998l) : ^Fronte ra et al.l 
j2000l') : JlSoderberg et al.l l l2004al'): ISakamoto et"all | |2004^ : ' Vandcrspek et _al.l J200^~3SLS^ 
1W04l : ^TM06l: ^this workp^Burrows et al.l l|2005l 'l: " 
1I2OO6I) : l ^Soderberg et al.l 1I2OO8I I. 



^ ^Sakamoto et al. I l|2008l 'l: A Campana et al.l 



Finally, we note that, given that SN 2008D manifests 
as an ordinary type Ibc SN, and that rate estimates for 
SN 2008D-type X-ray outbu rsts are consistent with the total 
rate of type Ibc supernovae ( Soderberg et ahlboosl ) , it seems 
likely that J5x,iso ~ 2 x 10*^ erg represents a lower bound on 
the prompt X-ray energy output of any GRB-SN. 



4 CONCLUSIONS 

The expanding X-ray halos observed by XMM-Newton in 
association with GRB 031203 provide a rare opportunity to 
study the prom pt soft X-ray properties of a nearby GRB- 
supernova (|VOJ). We have undertaken a detailed investiga- 
tion of these h alos in hopes of clarifying the nature of the 
"X-ray blast" (jwod ) whose dust-scattered X-rays are re- 
sponsible for the halos. 

In agreement with previous authors (|V04l : |W04|: Iwod : 
TMod ) . we find that the properties of the halos are consis- 
tent with an X-ray blast occurring simultaneously with the 
gamma-ray burst and subsequently scattering off dust sheets 
at distances d\ = 871.7 ± 9.7 pc and ^2 = 1388 ± 14 pc away 
in the plane of the Milky Way. In addition, we discover a 
third expanding halo that reveals scattering by a third dust 
sheet at distance ds = 9.94 ± 0.39 kpc; the presence of this 
third halo is a robust feature of our model as applied to the 
independent EPIC-MOS and EPIC-PN datasets, and due 
to the warping of the Milky Way disk, its location is well 
within the densest Hi regions of the outer disk. 

Our constraints on the timing of the X-ray blast, to = 
11 ± 417 s by comparison to the GRB 031203 trigger time 
(to < 545 s at 90%-confidence) , a re com parable in precision 
to the previous best constraint ( WOd l while suggesting a 
greater degree of synchronization between the GRB and the 
X-ray blast. Combined with upper limits from INTEGRAL 
observations of GRB 031203, which extend to t ~ 300 s af- 



ter trigger, our findings significantly restrict the parameter 
space of allowed models for the X-ray blast. 

We explore the implications of our findings, in reference 
to the properties of all other spectroscopically-confirmed 
GRB-SNe, and conclude that two alternative interpreta tions 
seem possible. On the one hand (jGhisellini et al.ll2006t ). the 
X-ray blast may be the signature of a high-energy "shock 
breakout" event, as observed from XRF 060218; in this case, 
the event should exhibit roughly x2 faster evolution and 
25% to 200% greate r X-r ay energy than the XRF 060218 
event. Alternatively ( WOSh . the X-ray blast may be the sig- 
nature of an X-ray fiare that occurred during the inter- 
val 300 s < to < 545 s; if so, the ratio of X-ray flare to 
gamma-ray burst energies for GRB 031203 would be high 
but not unprecedented. 

We point out that observations of XRF 060218 by Swift 
and IPN experiments allow for the presence of a "precur- 
sor" gamma-ray burst as luminous as (or more hkely, some- 
what less luminous than) GRB 031203. In this case, a shock- 
breakout origin for the GRB 031203 X-ray blast would make 
GRB 031203 and XRF 060218 "cosmic twin" explosions with 
nearly-identical high-energy properties. If GRB 031203 and 
XRF 060218 thus represent two facets of a single type of 
GRB-SN, detection and observation of future nearby GRB- 
SNe by Swift can be expected to yield additional examples 
of these events. 
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APPENDIX A: RING SURFACE BRIGHTNESS 

In this appendix we develop a semi-analytical treatment of the surface brightness of ring-type sources as observed with XMM. 
The results are incorporated into the numerical model of expanding dust-scattered halos which is used in the main text to 
analyse the soft X-ray properties of GRB 031203. 

We begin with the angularly-symmetric expression for the XMM point-spread function (PSF; Eq. [T]in the main text), 
expressed as a King function: 

sAr) = ^{l + r^/rl)-'', (Al) 

where r is the radial coordinate and the PSF parameters, core radius rc and scaling index a, are functions of the incident 
photon energy and exhibit slightly different functional forms for the MOS-1, MOS-2, and PN detectors; representative values 
for the parameters are Tc = 4.9" and a — 1.44. As expressed here, the surface brightness distribution is normalized to unit 
integral over < r < oo as long as a > 1. We note that the XMM PSF is known not to be angularly-symmetric and that this 
form is adopted strictly as a useful approximation to the actual PSF form. 

Figure lAll serves to define the coordinate system we use in our analysis. We seek to describe the surface brightness at 
the point x, due to a ring of radius R, as an integral over contributions from various angles, < (j) < 2n, where the angular 
coordinate is defined from the perspective of the ring centre. Owing to the angularly-symmetric nature of the PSF, we may 
take the point x to lie on the a;-axis, as shown, without loss of generality. 

The integral we wish to solve is thus expressed as: 

where the squared distance from the point x to the ring element at angle (f> can be derived by simple trigonometry as 
= R'^+x^-2Rx cos 0. Inserting this definition into the expression for the surface brightness, and factoring out all coefficients 
on the cos (p term, we find: 

sr(x) = ^ r,^("-i) {2Rx)-^ (^d+^±^ _ cos " d^, (A3) 
where the integral is now in a form that can be solved by use of one of the hypergeometric functions, 2F1. In particular, if we 
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Figure Al. (left) Coordinate definitions for analytical treatment of the surface brightness of ring-like sources, as observed with 
XMM. The ring has angular radius 7?, and we seek an expression for the surface brightness an angular distance x from the 
ring centre, which we can choose to lie on the a;-axis, as here, without loss of generality in an angularly-symmetric treatment. 
In analysing the contribution from a particular differential segment of the ring, at angular coordinate <f>, we define the distance 
from the segment position to the point x to be r. 



Figure A2. (right) Example surface brightness profile for a ring-like source, and comparison of our approximate approach 
to a full numerical integration. Upper panel: Surface brightness profile for a ring of radius R = 150" convolved with the 
azimuthally-symmetric approximation of the XMM point-spread function (Eg. IAl() . Two distinct approximations discussed 
in the text, Sr2{x) and Srn{x), are employed. The Sr2 approximation is used close to the ring centre (x < a;scp), and is plotted 
as a dashed green line beyond this domain. The Sm approximation is used for x > Xsop, and is plotted as a dashed red line 
near the ring centre; in the case shown, Xsup ~ 31.5". Lower panel: Residuals of the approximate approach, as compared to a 
full numerical integration of Eq. IA3I 



define A = (r^ + + x'^)/2Rx, then 

'(A-cos,^)-" d<^ = ^(A-l)-" ^F^C/2,a,l,^) +^(^+1)"" 2F1 ( V2, a, 1, ^), (A4) 



Sr2{x) = ^ r'J"-'^ (r'^+R' + x') 



1 + 



(A6) 





where we make use of the fact that yl > 1 in the present case (since r^, R, and x are all positive real numbers). The 
hypergeometric function 2-F1 has a series expansion that is valid and convergent over the full domain of interest; however, for 
the sake of computational speed we take two distinct approaches in our numerical calculations. 

In the first case, near the ring centre where x <^ R, we have 1, and the fourth argument of the 2F1 function in each 
term of the integral solution is small. In this case we expand each term to second order in the fourth argument: 

2^^1(72, a,l,z)^l + ^ + ^aia + l)z' + ■■■ (A5) 

so that, keeping terms to second order in 1/A throughout, the overall expression for the ring surface brightness becomes: 

a[a + 1)" 

This is our expression for the ring surface brightness in the second-order approximation. 

For x values near or beyond the ring, x ^ R, with R^ Vc still holding, we take a different approach. Writing x = R + A, 
we expand the integral expression for the ring surface brightness to second order in A/R: 

s™(^) = ^'-c^'^-'' (2ii^)-"£^ (l + ^ + ^"Cos0) d<l,. (A7) 

The integrand in this expression consists of a single, strictly positive, parenthetical term that is taken to the (—a) power. In 
the case of the XMM PSF model, a ~ 1.5 > 1. Thus, the integral will be dominated by contributions from the regime where 
the parenthetical term is near its minimum. Or, to take a more practical perspective, the bulk of the surface brightness will 
be contributed by the relatively nearby portions of the ring that have (j) ^ 1. Separately, we note that the integral limits may 
be changed from < </> < 27r to — tt < <j} < tt without efi^ect. 

We then make two further approximations, expanding to second-order in (j), coscf) ~ 1 — Va^^i s^ad increasing the now- 
symmetrical limits of integration from ±7r to ±00. The expression for the near-ring surface brightness then becomes: 

sUx) = ^ rl^-'^ {2Rx)- r f + ° #, (A8) 



where we have avoided expanding s: to i? + A in the {2Rx) term only. Making a change of variables to 1/' = (j)^ R^ /{A'^ + r-|) 
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leads to the simplified form: 



^rnix) = ^ r'J"-'^ {2Rx)-" ^/2 [^^^) J Jl + ^'V" d^, (A9) 

where the integral is now in a form that can be solved via the Euler gamma function: 

(i + ^r"d^-^^^#T^. (AlO) 

As a result, the final expression the ring surface brightness in our "near-ring" approximation is: 

. fli^S^l (l^) . (All) 



7r3/2y2 " ' ' V 2i?2 y V r(a) 

It remains for us to determine when to use each approximation. In our code we use one or the other approximation 
depending on the comparison of the value of x to the "separating value," a;scp: 



2ar ( ^l-(4a-l)(l + r2/r2) 



4a- 1 \ 2a 

with the test applying as follows: 



(A12) 



sr{x) = {'''^^''^' ^-^-p; (A13) 



An illustrative comparison of our immerical approach to a full numerical integration of Eq. IA3I is shown in Fig. IA2I The 
systematic underprediction of the true surface brightness in the Sm approximation is remedied when we renormalize the ring 
surface brightness over the active area of the detector, via a numerical integration. Overall, we find this level of accuracy to 
be sufficient to our purposes; moreover, as discussed in the main text, our use of a bootstrap Monte Carlo approach means 
that any inadequacies in our modeling should result in conservative estimates of parameter uncertainties. 
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